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Abstract 

We investigate a one dimensional chain of 2N harmonic oscillators in which neigh- 
boring sites have their energies redistributed randomly. The sites —N and N are 
in contact vi^ith thermal reservoirs at different temperature r_ and r+. Kipnis, 
Marchioro, and Presutti [18] proved that this model satisfies Fourier's law and that 
in the hydrodynamical scaling limit, when N oo, the stationary state has a lin- 
ear energy density profile 0{u), u & [—1,1]. We derive the large deviation function 
S{9{u)) for the probability of finding, in the stationary state, a profile 9{u) different 
from 9{u). The function S{9) has striking similarities to, but also large differences 
from, the corresponding one of the symmetric exclusion process. Like the latter it is 
nonlocal and satisfies a variational equation. Unlike the latter it is not convex and 
the Gaussian normal fluctuations are enhanced rather than suppressed compared 
to the local equilibrium state. We also briefly discuss more general model and find 
the features common in these two and other models whose S{9) is known. 



Keywords: Stationary non reversible states, Large deviations. Boundary driven 
stochastic systems. 
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1. Introduction 

The properties of systems maintained in stationary nonequilibrium states (SNS) 
by contacts with very large (formally infinite) thermal reservoirs in different equi- 
librium states are of great theoretical and practical importance. These arc argiiably 
the simplest examples of nonequilibrium systems to which the elegant, universal, 
and successful formalism of equilibrium statistical mechanics might hopefully be 
extended. A striking universal feature of equilibrium systems is the Boltzmann- 
Einstein relation according to which fluctuations in macroscopic observables, arising 
from the grainy microscopic structure of matter, can be described fully in terms 
of the macroscopic thermodynamic functions (entropy, free energy) without any 
recourse to the microscopic theory. In trying to develop a similar formalism for 
SNS we have to start with the fluctuations. There has therefore been much ef- 
fort devoted to developing a mathematically rigorous fluctuation theory for simple 
model SNS. This has led to some interesting recent results for conservative systems 
in contact with particle reservoirs at different chemical potentials [4-7, 10-12]. 

In particular it has been possible to obtain explicitly the large deviation func- 
tional (LDF) for some one dimensional lattice systems. The internal dynamics 
of these systems is governed by simple exclusion processes, symmetric (SEP) or 
asymmetric (ASEP), while the entrance and exit of particles at the two boundaries 
are prescribed by the chemical potentials, A±, of the right and left reservoirs. The 
LDF gives the logarithm of the probabilities of finding macroscopic density profiles 
p{u), where u is the macroscopic space variable, different from the typical values 
p{u); namely we have Prob(/9(u)) ~ exp{—NJ^{p)} where N is the number of lattice 
sites. 

In the symmetric case, the situation we shall be primarily concerned with here, 
the typical profile p{u) is given by the stationary solution of the diffusion equation 
dtp{t, u) = {l/2)du{Ddup{t, u)), u e [—1, 1] with boundary conditions p(±l) = p±- 
The values p± correspond to the densities in an equilibrium system with chemical 
potentials A±. The latter can be obtained by setting the chemical potential of both 
end reservoirs equal to each other, A+ = A_. We note that in this equilibrium case, 
the function J- is simply related to the free energy of the system. For p+ ^ p^ and 
constant diffusion coefficient D (that is density independent and spatially uniform) 
the profile p{u) is linear; this is the only case solved so far for the SEP. The results 
for the LDF of the SEP for this SNS contained some surprises. 

The most striking of these is non-locality: the probability of density profiles 
Pa{u) and Pb{u) in disjoint macroscopic regions A and B is not given by a product 
of the separate probabilities, i.e. the LDF is not additive. This is very different from 
the equilibrium case where the LDF is given (essentially) by an integral of the local 
free energy density for the specified profiles Pa{u) and Pb{u), and is thus automat- 
ically additive over macroscopic regions (even at critical points). Additivity is also 
true for the LDF of a system in full local thermal equilibrium (LTE), e.g. for the 
stationary nonequilibrium state of the zero range process. The microscopic origin 
of the non locality of the LDF for the open SEP lies in the 0{N~^) corrections to 
LTE which extend over distances of 0{N); N is number of lattice sites, which goes 
to infinity in the hydro dynamical scaling limit [1,24]. So while the deviations from 
LTE vanish in this limit their contributions to the LDF, which involves summations 
over regions of size N, does not 

The effect of these 0{N~^) corrections to LTE is already present at the level of 
Gaussian fluctuations about p{u). These were computed by Spohn in 1983 [24] who 
found that the contributions from the deviations from LTE made a finite contribu- 
tion to the variance of these Gaussian fluctuations, causing them to decrease, for the 
SNS of the SEP from their LTE values. The reduction in the variance of Gaussian 
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fluctuations can be recovered from the LDF by setting p{u) — p{u) + 7V~^/^0(u). 
In fact in [5, 11] it is shown that T{p) for the SEP dominates the LDF coming from 
the corresponding LTE state and therefore the fluctuations are suppressed. 

The above observations about the SEP raise many questions about the nature of 
the SNS of more reahstic systems. Do their LDF and Gaussian fluctuations behave 
similarly to those of the SEP? In particular, to what extent do the LDF for SNS 
play a "similar role" to free energies in equilibrium systems? In the absence of 
more solved examples it is difficult to answer these questions. It is therefore useful 
to flnd and investigate the SNS of other model systems for which the LDF can be 
found and compare them to that of the SEP. This is what we do in the present 
paper and then discuss the limited universality of the results. 

The SNS we consider here is a simple stochastic model of heat conduction in a 
crystal. It is well known, see e.g. [20,23], that harmonic chains do not obey Fourier's 
law of heat conduction. On the other hand, Kipnis, Marchioro, and Presutti [18] 
introduced a model of mechanically uncoupled harmonic oscillators in which near- 
est neighbor oscillators redistribute randomly their energy. This system is then 
coupled to thermal reservoirs at different temperatures and, thanks to the stochas- 
tic dynamics, the validity of Fourier's law is proven. In particular the stationary 
energy density 9[u) is a linear profile as in the SEP. We mention that a more 
sophisticated stochastic model of coupled harmonic oscillators has been recently 
investigated. The evolution is given by superimposing the Hamiltonian dynamics 
with a stochastic one in which two nearest oscillators randomly exchange momenta. 
This model has two conservations laws (energy and total length); the hydrodynamic 
limit is proven in [2] for the equilibrium case and in [3] for nonequilibrium, Gaussian 
fluctuations are analyzed in [16]. 

In this paper we consider the Kipnis-Marchioro-Presutti model, our main result 
is the derivation of the corresponding LDF, that we denote by S{6). It turns out 
that this function has both strong similarities and signiflcant differences from that 
of the SEP. Like for the SEP the LDF is nonlocal and yields Gaussian fluctuations 
about 9{u). Unlike the SEP, however, it is obtained by minimization, rather than 
maximization, of a "proto LDF" and the variance is increased compared to that 
obtained from LTE. Also in contrast to the SEP the LDF, S{6), is not convex. 
We discuss these similarities and differences in section |7| where we also give some 
generalization of our and previous results to a larger class of model systems. 

2. The model and main result 

Following [18] we consider a chain of one-dimensional harmonic oscillators lo- 
cated at sites x S [—N, A^] n Z =: Ajv and described by the canonical coordinates 
{qx,Px)- The oscillators are mechanically uncoupled so that the Hamiltonian of the 
chain is H = J2xi£An^p'x + Ix)/'^- The harmonic oscillators are however coupled 
by the following stochastic dynamics. Every pair of nearest neighbors sites waits 
an exponential time of rate one and then the corresponding oscillators exchange 
energy. More precisely, let {qy,Py), {qy+i,Py+i) be the canonical coordinates at the 
sites y, y + 1', when the exponential clock between y and y + 1 rings then the new 
values {qy,Py), {q'y+i,Py+i) are distributed according to the uniform distribution 
on the surface of constant energy 

Moreover the boundary site — A^, respectively -l-A^, waits an exponential time of rate 
one and then the corresponding oscillator assume an energy distributed according 
to a Gibbs distribution with temperature t_, respectively t+. All the exponential 
clocks involved in the dynamics are independent. 
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From a mathematical point of view it is sufficient to look only at the local energy 
given by the random variables S^x '■— {Px'^lx) /^i for which we get a closed evolution 
described by the following Markov process. The state space is Eat :— R^", an 
element of Sat is denoted by ^ := {^x , x G An}. The infinitesimal generator of the 
process is the sum of a bulk generator Lq plus two boundary generators i+ and 
L_ 

Ln ■.^N^[Lo + L^+L+] (2.1) 

in which we have speeded up the time by the factor N'^, this corresponds to the 
diffusive scaling. 

The bulk dynamics Lq is defined as 

AT-l 

Lq '■— Lx^x+1 
x=-N 

where 

Lx,x+if{0 f'dp [/(e(-'-+i)'^') - fiO] (2.2) 
Jo 

in which the configuration ^(^i^+i)'? is obtained from ^ by moving a fraction p of 
the total energy across the bond {x, a; + 1} to a; and a fraction 1 — p to a; + 1, i.e. 

{£.y if y 7^ a;, a; + 1 

Piix+(x+i) if y = x 

- P) iix + ^x+i) if y-a;+l 

The boundary generators L± are defined by a heat bath dynamics with respect 
to thermostats at temperatures r±, i.e. 

1 

i±/(0 := / dr — e-'-Z-i [fi^^^'n - /(O] 

in which the configuration ^i^'*" is obtained from by setting the energy at ±iV 
equal to r, i.e. 

lt:x,r \ \ if V ^ 

It — \ r \i y^x 

Note that we have set the Boltzmann constant equal to one. The process generated 
by lEH), denoted by ^(t), will be called the KMP process. 

We denote by w G [— 1, 1] the macroscopic space coordinates and introduce the 
space of energy profiles as M. :— {Q ^ Li([— 1, 1], dw) : Q{u) > 0}. We consider M 
equipped with the weak topology namely, 9n ^ iff for each continuous test func- 
tion (/) we have (j)) {0, </>), where (•, •) is the inner product in L2{[—1, 1], du). 
Given a microscopic configuration ^ G Sat, we introduce the empirical energy ttn{0 
by mapping ^ to the macroscopic profile 

JV 

[TTNmiu)--^ E ^-If^_^^+^l(") (2-3) 

I JV 2Af ' JV 2iV J 

X — — N 

note that ttn{0 G is a piecewise constant function. 

In the case when t_ = r_|_ = t it is easy to show that L^r is reversible with re- 
spect to the product of exponential distribution with parameter r, i.e. the invariant 
measure is given by the equilibrium Gibbs measure at temperature r, 

d^^NAO^ n ^e-«==/- (2.4) 

x=-N 
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When f e Ejv is distributed according to fiN,T then the empirical energy iTNiO 
concentrates, as TV — > cx3 on the constant profile r according to the following law 
of large numbers. For each S > and each continuous test function (j) — (j){u) 

lim /xjv..f|(7rw(C),0) - (t,<^)| >5)=Q (2.5) 

where r G is the constant function with that value. 

In this equilibrium case it is also easy to obtain the large deviation principle 
associated at the law of large numbers 1)2. 5|l . More precisely, the probability that 
the empirical energy t^n^^ is close to some profile £ different from r is 
exponentially small in N and given by a rate functional 5*0 

/^iv,r(^jv(C) -0) >=exp{ -A^5o(0)} (2.6) 

where ttn{0 ~ ^ means closeness in the weak topology of and x denotes loga- 
rithmic equivalence as N ^ oo. The functional 5*0 is given by 



50(6*) = / du 



1 



M-l-log^ 



= / dusoi9iu),0o) (2.7) 



where Oq = t is the constant energy density profile for t+ = t_ = r. The above 
functional can in fact be obtained as the Legendre transform of the pressure Ga(h) 

So{e) = sup [{e,h)-Go{h)] 

h 

where Go is defined as 

Go(/i) lim i-logE^„ fe^<''^-"(«») =- / du \og[l - t h{u)] (2.8) 

in which E^^ ^ denotes the expectation with respect to ^m,t- 

If T_ ^ r-|_ the process generated by L^r is no longer reversible and its invariant 
measure ij.n.t± is not explicitly known. Theorem 4.2 in [18] implies however the 
following law of large numbers. For each S > and each continuous (j) 

lim ^iN.r±(\{7TN{0,<t') -{0,(t>)\>s) ^0 (2.9) 

where 9 is the linear profile interpolating t_ and t_(-, i.e. 

-, , 1 — u 1 + u , , 

0{U)^T^^+T+^ (2.10) 

It is natural to look for the large deviations asymptotic for hn,t± ■ In the case of 
the symmetric simple exclusion process (SEP) this program has been carried out 
in [5,6,10,11]. The main result of this paper is an expression for the large deviation 
rate functional for fiN,T± analogous to the one for the SEP. The functional we obtain 
is nonlocal, as is the one for the SEP, but it turns out to be nonconvex while the 
one for SEP is convex. We mention that non convexity of the rate functional also 
occurs for the asymmetric exclusion process [12]. 

Without loss of generality we assume r_ < r+ and introduce the set 7^j_ := {t G 
Gi([-1, 1]) : t'(u) > , t(±1) = r±}, here t' is the derivative of r. Given 9 e M 
and r G 7^^ we introduce the trial functional 

Gi9, r) := fdum-1- log ^ - log ^^^] (2.11) 
J-i Lt(u) t{u) [T+-r_]/2J 

In this paper we show that the empirical energy for ^n.t± satisfies the large devia- 
tion principle with a nonlocal, nonconvex rate functional S{9) given by 

S{9)= inf g{9,T) (2.12) 
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that is we have 



Miv,r± {7tn{0 ~ ^) >; exp { - TV S{d) } 



(2.13) 



We note there is a very close similarity between (|2.12ll and the analogous result 
for the SEP, we emphasize however that in H2.12|l we minimize over the auxiliary 
profile r, while in SEP one needs to maximize. This is, of course, related to the 
non convexity of our S versus the convexity of the rate functional for SEP. It 
would be very interesting to understand this basic difference also in terms of the 
combinatorial methods in [10-12] besides the dynamical approach presented here. 

Given 6* G , wc show that the minimizer in H2.12|l is uniquely attained for some 
profile t{u) = t[9]{u); therefore S{9) = Q{6,t[6]). Moreover t[6'](u) is the unique 
strictly increasing solution of the boundary value problem 



We note that for 9 = 9 the solution of (|2.14|) is given by t[9] = 9 therefore 
S{9) = Q{9,9) = 0. On the other hand, by the convexity of the real functions 
M-l- 9 a; I— > a; — 1 — logx and M+ B x i-^ — log a;, for each 9 £ M and r S 75-^ we 
have Q{9, r) > hence S{9) > 0. By the same argument we also get that S{9) = 
if and only ii 9 = 9. This shows that the large deviation principle (|2.13|l implies 
the law of large numbers 1)2. 9|l and gives an exponential estimate as TV oo. We 
finally remark that the reversible case H2.7|l is recovered from (|2.11|I - H2.13|I in the 
limit T+ — r_ ^ which impose t(m) constant. 

Outline of the following sections. 

Our derivation of the rate functional S follows the dynamical/ variational ap- 
proach introduced in [4,5]. We look first, in Section |31 at the dynamical behavior 
in the diffusive scaling limit in a bounded time interval [0, T]. In particular, we ob- 
tain a dynamical large deviation principle which gives the exponential asymptotic 
for the event in which the empirical energy follows a prescribed space-time path. 

In Section 21 we introduce the quasi potential, it is defined by the minimal cost, 
as measured by the dynamical rate functional, to produce an energy fluctuation 9 
starting from the typical profile 9. By the arguments in [4,5], the quasi potential 
equals the rate functional S{9) of the invariant measure ^n,t±- A mathematical 
rigorous proof of this statement for the SEP is given in [7] . As discussed in [4, 5] , 
the quasi potential is the appropriate solution of a Hamilton-Jacobi equation which 
involves the transport coefficients of the macroscopic dynamics. The derivation of 
the functional S is then completed by showing that (|2.12|l is the appropriate solution 
of this Hamilton-Jacobi equation. As in the case of the SEP we are also able, by 
following this dynamical/ variational approach, to characterize the minimizer for the 
variational problem defining the quasi potential; this path is the one followed by the 
process, with probability going to one as iV ^ oo, in the spontaneous creation of the 
fluctuation 9. In Section^we also show that the functional S is not convex, obtain 
its expression for constant profiles 9, and derive an additivity principle analogous 
to the one for simple exclusion processes obtained in [11, 12]. 

In the remaining part of the paper we discuss some extensions of the previous 
results. In particular, in Section we discuss the KMP process in higher space 
dimension, d > 1, and obtain an upper bound for the quasi potential in terms of 
the local equilibrium one. We note that for the SEP it is possible to prove [5, 6] 
an analogous lower bound. We also discuss the Gaussian fluctuations around the 





(2.14) 



when 9 is kept fixed. 
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stationary profile 9] as for the SEP [5,9-11,24] the correction due to nonequihb- 
rium is given by the Green function of the Dirichlet Laplacian. In particular, this 
correction is non local; as in the case of the SEP, this is due to the long range 
correlations [18]. However, for the KMP process, the nonequilibrium enhances the 
Gaussian fluctuations while in the SEP it decreases them. As the covariance of the 
Gaussian fluctuations equals the inverse of the second derivative of S{6) at 9, the 
enhancement of Gaussian fluctuations corresponds to the upper bound of S{9) in 
terms of the local equilibrium functional. In the analysis in [18] a crucial role is 
played by a process, in duality with respect to the KMP process, in which the local 
variable at the site x takes integral values. In SectionElwe discuss briefly the large 
deviations properties of this dual model and obtain the expression for the large 
deviation functional. Finally in Section |7| we discuss the derivation of the large 
deviation functional for generic one-dimensional nonequilibrium symmetric models 
with a single conservation law. We obtain a simple condition, which is satisfied by 
the zero range process, the Ginzburg-Landau dynamics, the SEP, the KMP process 
and its dual, that allows the derivation of the large deviation function by means of 
a suitable trial functional. Even when this condition fails to hold, it yields a simple 
criterion to predict the enhancement/supression of the Gaussian fluctuation in the 
stationary nonequilibrium state with respect to the full local thermal equilibrium. 

The discussion in this paper will be kept at the physicists level of mathematical 
rigor. However, for the more mathematically inclined reader, we shall point out 
the main differences and technical difficulties with respect to the case of the SEP, 
which has been analyzed in full mathematical rigor [6]. 



3. Macroscopic dynamical behavior 

In this Section we consider the KMP process in a bounded time interval [0, T] un- 
der the diffusive scaling limit. We discuss the law of large numbers (hydrodynamic 
limit) and the associated dynamical large deviations principle for the empirical 
energy l|2.3(l . 

Given a continuous strictly positive energy profile 9 G C([— 1, 1];M+), we denote 
by i'^ the probability on Sjv corresponding to a local equilibrium distribution 
(LTE) with an energy profile given by 9. It is defined as 

N 
x=-N 

where 



Given two probability measures z^, on S^r we denote by h{i'\fj,) the relative entropy 
of I' with respect to /i, it is defined as 



We shall consider the KMP process with initial condition distributed according to 
the product measure for some energy profile 9q. A straightforward computation 
then shows there exists a constant C (depending on 6*0) such that for any N we 
have the relative entropy bound 

hii^lW^)<CN (3.1) 



where 9 is the stationary energy profile H2.10|l . By the weak law of large numbers 

< 



for independent variables we also have that i^^ is associated to the energy profile 
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lim 

N^oc 



0Q in the following sense. For each S > and each continuous cfi 

hm iy^J\{nN{0,<l>) - {Oo,^)\ >S)=0 (3.2) 

We remark that for the SEP it is possible (and convenient, see [6]) to consider 
deterministic initial conditions. For the KMP process, as the "single spin space" 
M-l- is not discrete, such initial conditions do not satisfy the entropy bound H3.1|l . 
which is required in the standard derivation, see e.g. [17,25], of the hydrodynamic 
limit. For this reason we have chosen the initial condition distributed according to 
the product measure I'g^. On the other hand, by the method developed in [21], it 
should be also possible to consider deterministic initial configurations. 

We denote by Pj^w the distribution of the KMP process when the initial condition 

is distributed according to Vg^^. The measure P^^n is a probability on the space 

Z?([0, r]; Eat) of right continuous with left limit paths from [0,r] to S^r. The 
expectation with respect to P^n is denoted by E^n . 

3.1. Hydrodynamic limit. Equation H3.2|) is the law of large number for the 
empirical energy at time t = 0; the hydrodynamic limit states that for each macro- 
scopic time t G [^iT] there exists an energy profile 9{t) such we have the same law 
of large numbers 

,«(|(7rAr(C(t)),(/.)-(0(t),(/.)| ><5) =0 (3.3) 

Furthermore, we can obtain the energy profile 9(t) by solving the hydrodynamic 
equation. For the KMP process (as for the SEP) this is simply the the linear heat 
equation with boundary conditions r±, i.e. 6{t) =^ 6{t,u) solves 

dt0{t) = \A0{t) 

9it,±l) ^ r± (3-4) 

where A is the Laplacian. Note that the stationary profile 6 in H2.10|l is the unique 
stationary solution of (|3.4|l . 

We give below a brief heuristic derivation, which is particularly simple for the 
KMP process, of the hydrodynamic limit. We refer to [14,15] for a rigorous proof in 
the case of the so called gradient nonequilibrium models with finite single spin state 
space; the extension to the KPM process should not present additional problems. 

Let (/) be a smooth function whose support is a subset of (—1, 1); from the general 
theory of Markov processes, we have that 

|e,« {{nN{m\ 0)) = V„ (LNinNim). 0)) (3.5) 

Since the support of is a strict subset of [—1,1], only N'^Lq contributes to 
L N {t^ N {£.it)) ■, <t>)- A simple computation shows that, when y ^ iiV, 

Lo^y = \[iy^i+iy+i-'^iy\ (3-6) 

we thus get 

^2 i.x/N+l/(2N) 

LN{TTNim),<f>) = ^ [Cx-i(i)+Cx+i(i)-2C.(i)] / du^{u) 

^ ™^»__ Jx/N-1/(2N) 



xeA„ Jx/N-l/i2N) 
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here An(I>{x/N) N^[(l){{x - 1)/N) + 0((a; + 1)/A^) - 2(j){x/N)] is the discrete 
Laplacian. The first step above comes from 13.6|l and H2.3|l . the second step from 
discrete integration by parts and last step from the regularity of (p. 

We have thus obtained the weak formulation of H3.4|l : it remains to show that 
also the boundary condition d{t, ±1) = t± is satisfied. For this we need to use the 
boundary generators N'^L±. These are Glauber like dynamics accelerated by iV^ 
so that the energy has well thermalized to its equilibrium value. We get 

(e±iv(i))«r± (3.7) 



A standard martingale computation shows that, with a negligible error as — + 
oo, 7rAr(f(t)) becomes non random. Wc can then remove the expectation value in 
the previous equations and get H3.3|l . 



3.2. Dynamic large deviations. We want next to obtain the large deviation 
principle associated to the law of large number (|3.3|l : more precisely we want to 
estimate the probability that the empirical energy 7rjv(^(i)) does not follow the 
solution of (|3.4(l but remains close to some prescribed path tt = 7r(t, u). This 
probability will be exponentially small in N and we look for the exponential rate. 
We follow the classic procedure in large deviation theory: we perturb the dynamics 
in such a way that the path tt becomes typical and compute the cost of such a 
perturbation. 

Let H = H{t, u) be a smooth function vanishing at the boundary, i.e. H{t, ±1) = 
0. We then consider the following time dependent perturbations of the generators 
Lx,x+i in (ES 

Jo 

Note that we have essentially just added a small drift N^-^VH{t, x/N) in the energy 
exchange across the bond {x,x + 1}. We denote by the distribution on the 

path space _D([0, Tl; Eat) of this perturbed KMP process. As before is the 

expectation with respect to . 

The first step to obtain the dynamic large deviations is to derive the hydrody- 
namic equation for the perturbed KMP process. We claim that for each t G [0,T], 
each continuous 4>, and each 5 > we have 



lim 



pf« (\{'^N{m).<t>) - m.<t>)\ >5)^Q (3.8) 

where Q(t) = 0{t^u) solves 




(3.9) 
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The argument to justify l|3.9() is similar to the previous one. Including the effect of 
the perturbation, the computation following H3.6|l now becomes (as before y ^ ztN) 



dp ^[H(t,y/N)-H(Uy+l)/N)m,+,~(l-p)^,] ^^^^^ ^ _ ^ 



[i/(f,(y-l)/7V)-F(t,y/7V)]- 



2 

+ [i? (i, y/N) - Hit, {y + l)/N)] + + ^g+i 

As before, we consider a smooth function (j) whose support is a strict subset of 
(—1,1); then only N'^Lq contributes to L^{TTN{S^{t)),(f)) and we get 

tHi fcf,\\ J,\ ^ 1 Ar2 ,/ /An f ix-ljt) + ga; + l(t) - 2^a,(t) 

^Ar(7rjv(U^)),W -772^^ ^(a;/^^)! ^ 

+ [H{t, (x - l)/N) - H{t, x/N)] 

+ [H{t, x/N) - H{t, {x + 1)/N)] -^^W^^+iW+^^W'+^^+iW' 

3 

+ — 2^ ^ VArff (t, a: /TV) VAr0(a;/7V) 

where VnHx/N) := N[f{{x + 1)/N) - f{x/N)] is the discrete gradient. In the 
above computations we just used Taylor expansions and discrete integrations by 
parts. With respect to the very simple case discussed before, we face now the main 
problem in establishing the hydrodynamic limit: the above equation is not closed in 
7rAr(^(i)), i.e. its right hand side is not a function of TrN{^{t)). In order to derive the 
hydrodynamic equation H3.9|l . we need to express —^x^x+i+^x+^x+i terms of the 
empirical energy 7rAr(^). This will be done by assuming a "local equilibrium" state, 
we refer to [6, 14, 15, 17,25] for a rigorous justification in the context of conservative 
interacting particle systems. 

Let us consider a microscopic site x which is far from the boundary and introduce 
a volume V, centered at x, which is very large in microscopic units, but still infini- 
tesimal at the macroscopic level. The time evolution in V is essentially given only 
by the bulk dynamics N'^L^; since the total amount of energy in V changes only 
via boundary effects and we are looking at what happens after OiN"^) microscopic 
time units, we expect that the system in V has relaxed to the micro-canonical state 
corresponding to the local empirical energy ■jT]s[{^{t)){x/N). To compute this state 
let us construct first the canonical measure in V with constant temperature r > 0, 
namely the product measure diyv,Ti£,) ■= Yixev '^^^ e~'>^/'^. Let now mvM be 
the associated micro-canonical measure with energy density 9, i.e. 



xev 



We introduce the function a (6) defined by 



aie) := hm E^,_, ( - ^.^.+1 + + ^'+1) (3.10) 
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where we recall that Emve denotes the expectation with respect to the probability 
mv,e- By the equivalence of ensemble we can compute a{9) also as 

a{e) = E,, , ( - +C+ C+i) = 3 0^ 

According to the previous discussion, the system in the volume V is well approx- 
imated by a micro-canonical state with energy density nN{^{t)){x/N). As it is 
shown by the standard proofs in hydrodynamic limits, see e.g. [17,25], we can thus 
replace, for N large, + ^^.(i)' + ^x+i W with 3[7rjv(eW)(a:/A^)]'. We 

then obtain 

^E,« « i(7r^(e(0), A0) + {7TNm)'^H,\/<P) (3.11) 

which is the weak formulation of (|3.9|l . The arguments to show that the boundary 
conditions 9{t, ±1) = t± are satisfied and to remove the expectation value are the 
same ones as in the derivation of 1)3. 4|l . 

Let TT = Tr{t,u), {t,u) G [0,T] x [—1,1] be a given path. We recall that our 
task is to estimate the probability that the empirical energy TTNi^{t)) is close to 
7r{t) (short for n{t,u)). We write this probability in terms of the perturbed KMP 
process, namely 

P.- i^Nim) - At), t e [0, T]) = E^„ j (3-12) 

Equation 1)3. 9|l tells us for which H the path tt becomes typical for the perturbed 
KMP process. We thus choose H{t, u) so that 

\ H{t,±l) ^ 

which is essentially a Poisson equation for H (recall that tt is fixed). With this 
choice we have, for N large, P-^ {T^N{S,{t)) ^ '^{t)) ~ 1 and to derive the dynamical 

large deviation principle we only need to compute the Radon-Nykodim derivative 
dP «/dPf„. 

We consider first the case of a deterministic initial configuration G '^n- In 
this case, by a standard computation in the theory of jump Markov processes, see 
e.g. [17, Appendix 1.7] or [5, Appendix A]), we have 

^(e)=exp{-iVJ[^^j(e,i?)} 

where 

J[o,T]itH) := (^w(e(r)),H(r)) - {nr,{^o),H{0)) - dt {7rr,m),dtH{t)) 

Jo 

N-l i-T i-l 
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By Taylor expansion we then get 

r-T 



/ E Ut)^NHit,x/N) 

■^0 x = -N+l 

n 

^_jv W ^[^(i, -1 + 1/^) - H{t, -1)] 

1 /-^ 

+ 2 7^ dUjvW^^[^r(t,l)-i?(i,l-l/A^)] 

'^^2]V ^ 3 [VNH{t,x/N)\ 

By the same argument given in the derivation of the perturbed hydrodynamic equa- 
tion we can replace -^xit)^x+i{t) + ^xit)^ + ^x+iit)^ by 3[Tr Ni^it)){x/N)]\ 
Recalling that in H3.12|l there is the indicator of the event in which 7r7v(^(t)) is close 
to 7r(t), we get 

J[o,T]{^,H) « J[o.T]W - (^(r),i/(r))-(^(0),il(0))- / dt{7r{t),dtHit)) 

Jo 

- \ j\t {nit), AHit)) - i £dt (7^(^)^ [VHit)]') 

dtT^VH{t,~l) + ^J^ dtT+VH{t,l) 

where we used the fact that the value of tt is fixed at the boundary, 7r(i, ±1) = 
T±. Recalling that the perturbation H has been chosen as the solution of H3.13|l . 
integration by parts shows that 

J[o,T](vr) dt {VH{t), n{t)^S/H{t)) (3.14) 

To complete the derivation of the dynamical large deviation functional, we only 
need to consider the fluctuations of the initial condition. Recalling that we have 
chosen the initial condition distributed according to the product measure , a 
straightforward computation on product measures (the one carried out in H2.6|l - 
(EiHI)) shows that 

<('riv(e) 7r(0)) X exp{ - iV5o(^(0)|eo)} 

where S'o(7r(0)|^o)7 which represents the contribution to the dynamic large deviation 
from the initial condition, is given by 



^0(^(0)1^0) = / du 



1 



1 



7r(0, m) 7r(0,u 



1 - log ■ 



%iu) " 0oiu) 



(3.15) 



By collecting all the computations performed we finally get the dynamical large 
deviation principle 

P,«(7rAr(e(t))~^(t),ie [0,r]) xcxp{-7V/[o.T](7r|0o)} (3.16) 

where 

I[o,T]{n\9o) = Soiirmeo) + J[o,T]M (3.17) 
We note again that S'o(7r(0)|6'o) represents the cost to create a fluctuation at time 
zero whereas J[o,t]{'^) represents the dynamical cost to follow the path 7r(i) in 
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the time interval [0, T]. In the case of deterministic initial conditions, as the one 
discussed in [6] for the SEP, we would have S'o(7r(0)|6'o) — +oo unless 7r(0) = 6o- 

3.3. Remarks. We conclude this section with some remarks on the rigorous deriva- 
tion of the dynamical large deviation principle H3.16|l . The probability estimates 
needed are (not surprisingly) more subtle than discussed here. In fact, while in 
the proof of the hydrodynamic limit it is enough to show that we can replace 
-^x{t)^x+i{t) + Q{t) + il+^{t) by 2,['KN{i{t)){x/N)f with an error vanishing as 
N oo, in the proof of the large deviations we need such an error to be of o{e^'~''^). 
This can be proven by the so called super exponential estimate, see [17, 19], which 
is the key point in the rigorous approach. This estimate has been extended to the 
non equilibrium SEP in [6]. For the KMP process there is the additional com- 
plication of a unbounded single spin space. In [13] the dynamical large deviation 
principle is proven for the Ginzburg-Landau model; however for the KMP process 
the situation is more troublesome because the mobility tt^ is unbounded and the 
reference measure is only exponentially decaying for large tt. There is also another 
technical point which requires some care. In the usual proofs of large deviations 
from hydrodynamic behavior, one first obtains the lower bound for a neighborhood 
of strictly positive smooth paths tt and then uses approximations arguments to 
extend the lower bound to any open set. The approximations arguments used for 
the SEP, see [17, 19] for the equilibrium case and [6] for nonequilibrium, take full 
advantage of the fact (special for the SEP) that J[o.t] (""j is a convex functional. In 
order to prove the dynamic large deviation principle for the KMP process a more 
robust approximation method, possibly analogous to the one in [22], is required. 

4. The quasi potential and its properties 

In this Section we introduce the quasi potential, which measures the minimal 
cost to produce a fluctuation of the energy profile in the stationary state, and shows 
that it can be obtained by solving the one-dimensional non linear boundary value 
problem (|2.14|l . We also characterize the most probable path followed by the KMP 
process in the spontaneous creation of such a fluctuation. We finally show that the 
functional S is not convex and derive an additivity principle analogous to the one 
in [11,12]. 

Given T > and a strictly positive smooth 6 € A4, we introduce the set of 
energy paths which connect 9 to 6 in a time interval [— T, 0], i.e. we define 

Sg j. {tt = TT{t, u) : n{-T, u) = e{u) , 7r(0, u) = e{u)} (4.1) 

where we recall that the stationary energy profile 9 has been defined in H2.1()|l . 
Paths TT e Se^T must also satisfy the boundary condition 7r(i, ±1) = t±; in fact it 
can be shown [6] that J[-t,o] (^) = +oo if the path tt does not satisfy this boundary 
condition. The quasi potential is then defined as 

V{9):^M inf J[_T,o](^) (4.2) 

where we recall that the functional J is given in H3.14() . By the general arguments 
in [5] , see also the rigorous proof in [7] for the SEP, we have that the rate functional 
S{9) for the invariant measure ^j,n,t±, see ()2.13|l . coincides with the quasi potential, 
i.e. S{9) = V{9). 

4.1. Solution of the Hamilton— J acobi equation. Recalling that the perturba- 
tion H in (|3.14|) solves (|3.13|) . the variational problem 1)4. 2|l consists in minimizing 

the action that correspond to the Lagrangian 

C{9,dt9) = l{V-'{dt9 ~ ^A9), ^V-\dt9 ~ ^A9)) (4.3) 
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The associated Hamiltonian is 

n{e,H) sup {(iJ,C) -mC)} = \{^H,e^\JH) + \{H,M) (4.4) 

Noting that V is normahzed so that V{0) = 0, we obtain, by a classical result in 
analytic mechanics, that ViO) solves the Hamilton-Jacobi equation H(0, ^) = 0, 
i.e. 



5V ^ 5V\^ iSV 







(4.5) 



where 5V /56 vanishes at the boundary and 0{±.l) = t±. We look for a solution of 
(|4.5|) in the form 

for some function r — t[9]{u) to be determined satisfying the boundary conditions 
r(±l) = T±. By plugging H4.6|l into H4.5|l and elementary computations, analogous 
to the ones for the SEP discussed in [5] , we get 



'vi 

r 



-(v(.-r),5:) + (<^,(«^-r^)) (4.7, 



= 



T 

Therefore a solution of (|4.7|l is obtained when t satisfies the non linear bound- 
ary value problem (|2.14() . Let us denote by t[6] the solution of (|2.14|l : recall the 
definition H2.11|l of the functional G{d, t) and that, since H2.14|l is the associated 
Euler-Lagrange equation for fixed 6, we have [SQ/St] {9,t[9]) =0. By a direct 
computation we then get 



(4.8) 



which shows that V{9) — Q{9,t[9]) is a solution of the Hamilton-Jacobi equation 
H4.5|l . To complete the derivation of H2.13(l we next show that V{9) meets the 
criterion in [5, §2.6], i.e. it is the "right solution" of the Hamilton-Jacobi equation, 
and that the infimum in H2.12|l is uniquely attained for t = t[9], the solution of 
(EH- 

4.2. The exit path. The characterization of the optimal path for the variational 
problem (|4.2|l can be carried out according to the general scheme in [6]. Let V{9) = 
Q{9,t[9]) and 7r(i), t G [— T", 0] a strictly positive smooth path such that 7r(0) = 9. 
By using that V{9) solves Hamilton-Jacobi equation (|4.5(l . a simple computation 
shows that 



J[_T.o](7r) = 2 



dt ( V 



-T 



9,^-iA^ + v(7r2v^(^)) -v(7r2V^(7r 



9t7r- ^A7r + VItt^V — (tt) 



V{9) - V{n{-T)) + 



59 



5V_ 
19 



dt ( V 



-T 



dtn - ^Att + v(^7r^V^(7r) 



-A. + v(.^V-(.)^ 
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Since the last term above is positive, the optimal path tt* for the variational problem 
(|^ solves 

d,7:* = Un* v(in*fvZin*)) = -^An* + v(-^^Vr[n*]) (4.9) 



2 V ' Se' 'J 2 MT[7r*])2 

where T[7r*] = T[TT*]{t,u) denotes the solution of (|2.14() with 6 replaced by 7r*(t). 

Let us denote by 9*{t) = 7r*(— i), t G [0,T], the time reversed of the optimal 
path TT*. It is then not difficult to show, see [5, Appendix B] for the analogous 
computation in the case of the SEP, that 9* (t) can be constructed by the following 
procedure. Given 9 = 7r*(0) = 9*{0), first let tq = t[9] be the solution of H2.14|l . 
then solve the heat equation with initial condition tq, i.e. let T{t) be the solution of 

tq{u) 
and finally set 

Since r(t) ^ 6* as t ^ c» we get 7r*(-T) ^ 6* as T ^ oo, hence V{'!t*{-T)) 
V{9) — 0. The identification of the solution of the Hamilton-Jacobi equation with 
the quasi potential follows from the characterization of the minimizer tt* obtained 
before. In particular V{9) satisfies the criterion discussed in [5, §2.6]. 




4.3. Solution of equation (|2.14|) . The existence of a solution for the nonlinear 
boundary value problem H2.14|l can be proven by the same strategy as in [6, 11]. 
We write H2.14|l in the integral-differential form 

r r [Tiw)-9iw)]T'iw) ^ 
av exp < / aw -. — > 

1 U-i nwy J 



t{u) = t_ + (r+ - r_) —J- 



..exp^ / ,JrH-OH]r'M 
r(u))2 

iferen 

t{w) ~ 9{'w)]t'{w) 



-1 

Then a solution of 12.14|l is a fixed point of the integral-differential operator ICg [t] 
defined as 



ICeir] (w) := t_ + (r+ - t_) 



1 



dv exp < / dw 



r(w)2 



..exp^ / ,JrH-OH]r'H' 
t{wY 



1-1 

We consider the case in which 9 is bounded, namely we assume that \\9\\ := 
sup„|0(w)| < +00. Recalling that t must be strictly increasing and such that 
t(±1) = T±, we get 

\\9\\t' ^ {r-9y ^ r' 



which yields 

It is now easy to show, see [6] for more details, that for each G A^, the operator 
/Ce[T] maps a compact convex subset of 7^^ into itself. Hence, by Schauder's fixed 
point theorem, we conclude the proof of the existence of solution to H2.14|l . 

Uniqueness of solution to H2.14|l can also be proved with a slight variation of the 
argument in [11]. Let us consider two different increasing solutions of (|2.14() ti(u) 
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and T2{u). If t{(— 1) = T2(— 1) then uniqueness of the Cauchy problem impHes 
Ti = T2. On the other hand, if r{(— 1) > T2(— 1) > then we denote by u the 
leftmost point in (—1,1] such that ti(u) = T2(u). The point u exists because 
Ti(l) — T2(1) moreover we have that t[{u) < T2(u). From H2.14|l we get 



d z T 
du \t' 



which integrated gives 



t'{u) ^ r'(-l) 



we then deduce 
ri(u) 



T2(m) 



Tiiu) r^{u) T[{-1) T^(-l) 



1 ^(l') 



e{v) 



1 



-Ti(w) r2(w)- 

Since r{(— 1) > T2(— 1) and ti(w) > T2{v) for v G [— 1,m], the right hand side above 
is strictly negative. Recalling that ti{u) = T2(u) we get t[{u) > T2(m), the desired 
contradiction. 

In order to prove that the infimum in (|2.12|l is uniquely attained for t = t[0], 
the solution of H2.14|l . we perform the change of variable t = e'^. We then get the 
functional 

"1 ^ V'{u) 



g{e,ip) ■=g{d,e'^) = / du e{u)e-'^^'''> -i-\oge{u) ~\og- 



1 ^ [r+ 
which is strictly convex in ip; this trivially implies the claim. 



r_]/2 



(4.10) 



4.4. Non convexity of the quasi potential. As we mentioned before, in the case 
of the SEP the quasi potential can be obtained by a variational problem analogous 
to (|2.11|l where one maximizes over the auxiliary profile [5,6,10,11]. In such a case, 
since the functional Q{6,t) is convex in 9 for fixed r, the rate functional S{6) is 
trivially convex in 0. In the case of the KMP process we need instead to minimize 
over the auxiliary profile r, therefore there is no reason to expect S{9) to be convex. 
We now show, by an explicit computation, that the rate functional is indeed not 
convex. 

We mention that non convexity of the rate functional S has been shown for the 
asymmetric exclusion process [12]; in that case however the functional is degenerate, 
in the sense that there are infinitely many profiles for which S vanishes. Therefore 
the mechanism of the non convexity is somehow different from the one in the KMP 
process, where S{d) vanishes only at 9. 

To prove the non convexity of the rate functional S we shall exhibit profiles 9 
and g so that, by choosing e small enough, we have 



S{9)^S[^[9- 



£9] + l[0-£9] 



>lsio 



eg) + ls{9 



eg) 



(4.11) 



Let T = t[9\ be the strictly increasing solution of the boundary value problem 
()2.14|l . then by using (|2.12|l for e small enough and any profile / vanishing at the 
boundaries, /(±1) — 0, we have 



S{9 + eg) < 



0(9 + eg, T 
g{0,r) 



/>{(^o(f-^)-^} 



du 



.,_,-{( 

where we brutally Taylor expanded H2.11() 



9 



.9f , if? 



} + o{e' 
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Since S{6) — Q{6,t), the inequality (|4.11() will follow if we show that, for an 
appropriate choice of / (recall that t = t[9] is the solution of H2.14|l ) we can make 
the quadratic term in the previous equation strictly negative, i.e. 

/>{(4-i){;4-f^-i{|}<° 

Let us introduce the function 



h{u) := 



+ ^(r+-r_)(l-M6) if -l<y<-l/2 
+ ±{r+ - T^){u - 1) if -l/2<u<l 



0{u) ^ Hu) 1 - hiu)jj^ (4.13) 



Note that h E C^([— 1, 1]), h{±l) = t± and h is strictly increasing. We choose the 
profile 9 as 

h"{u) 

'h'{uy 

Note that 6* > 0, i.e. 6 is an allowed profile in A^; the corresponding solution of the 
boundary value problem (|2.14l) is t[0] = h. We further choose f{u) = (1 — v?')h'{u) 
(note that / vanishes at the boundaries as required) and g = fO'^/h'^. With the 
above choices the left hand side of (|4.12|l equals 

f2 f2n2 (fl\2. 



J jii^ - u(l - u2)^| 



< 

This completes the proof of (|4.12|) and therefore of the non convexity of the rate 
functional S. 

4.5. The rate functional on constant profiles. Here we show that for constant 
profiles the boundary value problem (|2.14|l can be integrated; the corresponding 
value of the rate functional S{9) can be expressed in terms of special functions. 

We use the variable (p — logr; we then have S{9) = Q{d,Lp[9]), where the func- 
tional G has been defined in H4.1()|l and (p[9] is the unique strictly increasing solution 
of the boundary value problem 

if" 



((^')2 (4.14) 

ip{±l) = l0gT± 

If we restrict to constant profiles 9 this equation can be integrated obtaining 

log^'[6i](M) = log</7'[e'](-l) +6l(e-'^[^l(") - —1 (4.15) 
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and from this 
From equation H4.15|l we obtain also 

and finally with a change of variables 



T+ - T_ 



= 2^ + — -log0-log(^'[0](-l)+log 2 



1 S- 



log 



s{e) ^ 2 



log 



dy 



dip e 



l + log(r+ -T_) 



(4.16) 



In particular for 9 large, from H4.16|l we deduce the asymptotic expansion 

S,., = 2{A,(,„,^_1),^},0(1) («7) 

Recall that the equilibrium functional 5*0 is given in l|2.7|) and note that for 
constant and large values of the profile 9 we have So{9) « 29 /t. By comparing this 
with the expansion (|4.17() . we see that the leading order is the same but only the 
warmer thermostat matters, as it is quite reasonable from a physical point of view. 

As we showed earlier, the rate functional S is not convex. The restriction of S 
to constant profiles obtained in H4.16|l might however be convex; we do not have 
an analytic proof of the convexity of (|4.1t)|) , but rough numerical evidences suggest 
this is the case. 



4.6. An additivity principle. In [11] the rate functional S was derived for the 
SEP by combinatorial techniques. It was then shown that S satisfies a suitable 
additivity principle which allows to construct the rate functional for a system in a 
macroscopic interval [a, b] from the rate functional of subsystems in the intervals 
[a, c] and [c, b], here a < c < b. More precisely, in [11] is introduced a modified rate 
functional S[a,b] (''a, T),; 9) where Tq, ti, are the density at the endpoints and 9 = 9{u) 
is the density profile in [a, b]. The additivity principle is then formulated as 

S[a,b]{Ta,n;9) = sup {S{a,c]{Ta,Tc;9\^a,c\) + '5'[c,fc] (^c, Tf,; 6* [[^,6] ) } (4-18) 

where 9\]^a,c\i respectively 0[[c,6], denotes the restriction of the profile 0, which 
is defined on the interval [a, 5], to the subinterval [a,c], respectively [c, 6]. The 
additivity principle H4.18|l plays a crucial role in the derivation of the rate functional 
for the asymmetric exclusion process. In [12] the expression of the rate functional 
from this principle is then deduced. 

Here we show that the rate functional S for the KMP process satisfies an ad- 
ditivity principle analogous to (|4.18() . Here, however, we need to minimize on the 
midpoint value r^. This is due to the fact that in H2.12|l we need to minimize 
over the auxiliary profile r; a direct derivation of the additivity formula, as was 
done in [12] for the asymmetric exclusion process, would clarify the basic physical 
difference between the KMP process and the SEP. 

Let us consider the KMP process on the macroscopic interval [a, 6], here we 
denote the temperatures of the heat baths at the boundary by Ta, Tf,. We then let 
S[a,b]{Ta,Th]9) be the corresponding rate functional and introduce 

S[a.b] (Ta, n;9)^ S^aM [Ta, u; 9) - (b ~ a) log (4.19) 
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by using (|2.11() and (|2.12() we then get 



J a 



■b 



9{u) 



inf 



r : 



du 



It{u) 



1 - log 



t{u) 



logT'(u) (4.20) 



r(a)=Ta, T(b)=Tt, 



where the infimum is over the strictly monotone auxiliary profiles t(m), u € [a.,b]. 
We then get the additivity principle for the KMP process: 



It is not difficult to show, see [11], that the expression (|2.11(l - (|2.12|l for the rate 
functional follows from the additivity rule (|4.2HI . 



4.7. Remarks. We again conclude with a few mathematical remarks. We have dis- 
cussed existence and uniqueness of (12.14(1 only for bounded profiles 9; the extension 
to 9 € M should be however straightforward. Let V{9) be the quasi potential as de- 
fined by the variational problem 1(4. 2|l . Since the optimal path tt* has been explicitly 
constructed, the rigorous proof of the upper bound V{9) < infi-gr^^ Q{9, r) should 
be carried out as in [6]. The proof of the lower bound V{9) > 'mireT^^ G{9,t) is 
instead more troublesome. The computations presented here essentially prove this 
bound for strictly positive smooth paths tt, but the argument in [6] to extend it 
to arbitrary paths in Sg^T takes advantages of the convexity (special for the SEP) 
of the dynamical rate functional J. This problem is of course related to the proof 
of the lower bound for the dynamical large deviation principle for any open set 
mentioned at the end of Section |21 The identification of the quasi potential V with 
the rate functional for the invariant measure S has been proven for the SEP in [7] ; 
although the strategy is of wider applicability, the technical points might require 
some extra effort. 

5. Higher space dimensions 

The KMP process introduced in Section|2can be easily generalized to the case of 
space dimensions d > 1. Let A be a smooth domain in R"^ and set A^r := Z'^ n A'^A. 
We then define the process on S^r := R^"as follows: every pair of nearest neighbors 
oscillators exchanges energy according to the rule described in Section [21 and every 
oscillator at a boundary site x is in contact with a thermostat at temperature 
f{x/N) for a fixed function f. 

Several computations of this paper can be repeated step by step when the model 
is not one-dimensional. In particular the hydrodynamic equation has still the same 
form ((3.41) with the boundary condition 9\gA= f and the dynamic large deviation 
functional J has the same form as ((3.14(1 . Formula ((4.21) as well as the Hamilton- 
Jacobi equation 1(4. 5() for the quasi potential holds in any dimension; we can still 
perform the change of variables ((4.6(1 and reduce 1(4.5(1 to 1(4.7(1 . However the solution 
of the boundary value 1(2. 14() does not give the quasi potential because, with this 
choice, the right hand side of ((4.6(1 is a functional derivative only if d = 1 . 

However, by analyzing the variational problem 1(4.2(1 . we derive an upper bound 
for quasi potential V{9) that holds in any space dimension. We also discuss here 
the Gaussian fiuctuations of the empirical energy when ^ is distributed according to 
the invariant measure. We shall obtain the covariance of the Gaussian fluctuations 
by expanding the large deviations functional S{9) around the stationary profile 9. 
We note that in the one-dimcnsional case the arguments are easier thanks to the 
more explicit form of S. 



S[afi]iTa,n;9) 




'.Ta,Tc;9\[a^c]) + S[c^b]{Tc,Tb;9\[r^^b])} (4-21) 
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5.1. Upper bound for the quasi potential. Let us denote by 9{u), u € A the 
stationary solution of (|3.4() with boundary condition 0{u) = t{u), u € dA. Note 
that, for generic boundary conditions f , the profile 9 does not have the simple form 
l|2.1()(l . Of course 9 is still the most likely profile for the empirical energy under the 
invariant measure. We also introduce the local equilibrium large deviation function 



J A 



du 



9iu)_^_^^^9iu) 



(5.1) 



_9{u) 0{u)_ 

and note that it coincides with the function defined in (|3.15() and it is thus the rate 
functional for the product measure I'g' . 

When ci = 1 we can use (|2.11ll and easily obtain the upper bound 



S{9) 



inf 



g{9,T)<g{9,9)^s,^{9) 



(5.2) 



For d > 1 we use a different strategy. Given a path tt = Tr(t, u) satisfying the 
boundary condition Tr{t,u) = f(u), u G 9A, a straightforward computation shows 
that 



J[-t,o]{t^) 



-T 




dt ( 









1 A 1 


dtTT Att 




dtTT - 


-Att 


L 2 J 






2 J 



dt(v 



Att 



dtTT 



1 



Att + VlTT^V 



,1 



Att - V 



5eq(7r(0)) - SM-T)) - 



1 
2 



dt 



-T 



dtTT + -Att + VItt^V- 
2 \ 9 

dtTT + ^Att + vfTT^vi 



The quasi potential is defined by the variational problem 14.2|l . Hence, to obtain 
an upper bound for V{9) it is enough to exhibit a path tt which connects 9 to 9. 
We choose 7r(t) = 9{—t) where 9{t) solves 



\M{t)+V 



(9{tfvl 
uedA 



dt9{t) 

9{t,u) = f{u) 
[ 9{0,u) = 9{u) 

We note that the path tt connects 9 to 9 since 9{t) ^ 9 as t ^ oo. By using the 
path TT in the above expression for J[_t,o] ("') ^-nd letting T — > oo we get 



V{9) < 5eq(0) - 



1 



" dt(^^,[TTit)-9r)<S,^i9) 



which shows that the upper bound V{9) < Scq{9) holds in any space dimension. 
We also note that the above inequality is strict unless 9 = 9. 

5.2. Gaussian fluctuations. In the case d = 1 we can follow step by step the 
argument given in [11] for the SEP. We consider a small perturbation, 9 = 9 + eQ, 
of the stationary profile 9, and consequently have t[9] = 9 + eT, where, to first 
order in e, H2.14|l gives 



{V9y 



-AT -T = -e 



(5.3) 
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The functional S{9) has a minimum at 9 so that its expansion in e is 

S{e) = S{e) + ^e^{0, C-^Q) + o (e') (5.4) 

The operator C is the covariance for the Gaussian fluctuations of the empirical 
energy under the invariant measure fiN,T±- Since S{9) ~ Q{9,t[9]), we get 

1 r[r(u)-e(^)p , [VT(u)]2- 



= jjn — m m^^^\ / 

where we used Taylor expansions, integrations by parts, T(±l) ~ 0, and (|5.3|l . The 
operator W is defined as 

From equation H5.3|l we get Q — —WAT and this implies 

/ W \ / W^^ 

^^'^"^^>-(^^'(W^'^)^<'''(W^ 

so that the covariance C is given by 

C ^{V6fW ^~9'^\^{SI9f {-AY^ (5.6) 

The first term above is simply the covariance of the Gaussian fluctuations of the 
empirical energy for local equilibrium product measure i/^, while the second term 
represents the contribution to the covariance due to the long range correlations 
in the stationary nonequilibrium state. As in the case of the SEP [5,9-11,24] 
this correction is given by (— A)~^, the Green function of the Dirichlet Laplacian. 
Since (— A)^^ > 0, for the KMP process this correction enhances the Gaussian 
fluctuations, while in the case of SEP it decreases them. We also note that, by 
exploiting the duality of the KMP process with the process we shall introduce in 
Section |H| the expression H5.6() for the covariance could be rigorously deduced as in 
the case of the SEP [9,24]. 

To obtain the covariance of Gaussian fluctuations in the case d > 1 we instead 
argue as in [5] . Let us introduce the "pressure" as the Legendre transform of the 
rate functional 5(0), i.e. 

G(h):^^M^{{9,h)-S(9)) (5.7) 
e 

We then get that G(/i) satisfies the Hamilton-Jacobi equation dual to 1)4.5(1 . i.e. 



where /i(m) satisfies the boundary conditions ft.(u)|aA = 0. 

Let us denote by Goq the pressure associated via ((5.7(1 to the local equilibrium 
functional S'eq; we look for a solution of 1(5.8(1 in the form 

G{h) = Geq(/i) + (g, h) + S h) + o(h^) (5.9) 

for some function g = g(u) and some linear operator B. From ((5.4(1 we get 

G(/i) = + i (/i,G/i) +o(/i2) (5.10) 
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hence the second derivative of G at /i = is the covariance of the density fluctua- 
tions. By comparing (|5.9|l to (j5.10(l we find 



C 



S^G, 



eq 



Sh^ 



+ B = 9'^l + B (5.11) 

/i=0 



By plugging H5.9|l into H5.8(l and expanding up to second order in ft,, it is not dilRcult 
to show, see [5] for the case of the SEP, that g — and 

{h,ABh)^-{h,\V9\^h) (5.12) 

The operator B therefore satisfies 

i[AB + BA] = -|V^|2 (5.13) 

See [24] for another derivation of this equation based on fluctuating hydrodynamic 
instead of large deviations. 

From H5.13|l we see that if is constant (this condition can be realized by a 
suitable choice of the thermostat f), the operator B has the kernel 

B{u, v) = |V0p(-A)-i(u, v) (5.14) 

where (— A)^^(u,w) is the Green function of the Dirichlet Laplacian in A. The 
interpretation of H5.11|l and (|5.14(l is as in the one dimensional case (|5.6|l : we note 
the fact that i? is a positive operator can also be obtained as a consequence of the 
bound S{9) < S'eq(6'). 

6. The dual process 

The analysis in [18] is based on a duality relationship between the KMP process 
and another process we discuss next in the one-dimensional case. The state space 
is Sat := N'^" where N := {0, 1, . . . } is the set of natural numbers. If C = {£.x , x S 
Aat} e Ejv, the value at the site x can therefore be interpreted as the number 
of particles at x. As for the KMP process, at each bond {x,x + 1} there is an 
exponential clock of rate one; when it rings the total number of particles + ^x+i 
is redistributed uniformly across the bond {x,x + 1}. Moreover the boundary 
sites ±iV evolve according to a heat bath dynamics with respect to a geometric 
distribution with parameter p±. More formally, the infinitesimal generator has still 
the form l|2.1() but now the bulk dynamics is defined by 



^x,x+if (0 . ^/ ^ [/(?(^'^+')''=) - /(O] (6.1) 

where the configuration ^C^'^+i)'*^ ig defined as 

{^v if y^x,x + l 

k if y = x 

Cx + ^x+1 - k if y = x + 1 

The boundary part of the generator is defined as follows 

oo 

i±/(0 E?'±(l -^^±)'[/(^^'^'') - /(^)] (6-2) 

fc=0 

where p± e (0, 1) arc the parameters of the reservoirs and the configuration ^^''^ is 
defined as 

^^xM I if y ^ ^ (6.3) 
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If = p- = p the dynamics is reversible with respect to the product of geo- 
metric distributions of parameter p, i.e. the invariant measure is 

t^NAo= n (6-4) 

By a computation analogous to (|2.6|l - (|2.8(l . it is easy to show that when ^ is dis- 
tributed according to /ijv.p then the empirical density 7rAr(^), which is defined as in 
(|2.3|) . satisfies a large deviation principle with the rate functional 

Soie)^ fdu{e{u)iog^M + [i + eiu)]iogj^^^] (6.5) 

where the parameter 9 is related to p by the relation 9 ~ J^kLo^P^^ ^ P)*^ ~ 
{l-p)/p. 

When ^ p_ the model is no longer reversible and the invariant measure /i n,p± 
is not explicitly known. In the sequel we shall assume p- > p+. We can repeat the 
computations done for the KMP process and get the hydrodynamic equation. This 
is still the linear heat equation with the appropriate boundary conditions, i.e. 

9{t,±l) = 9^ = (6-6) 
P± 

[ 9{0,u) = 9o{u) 

As before the most likely density profile 9 is the stationary solution of 1)6. 6|l . 

To obtain the dynamic large deviation principle we introduce a smooth function 
H = H(t, u) vanishing at the boundary and consider the following time dependent 
perturbation of the generator L^^x+i in l|6.1|l 

ix + ix+i + 1 J^^g 

X [/(e^"^"+^^''^) - /(e)] 



The hydrodynamic equation associated to this perturbed dynamics is given by 

dt9{t) = ^A9it) - v(0(t)[l + e(t)]Vff(t)) (6.7) 

with the same boundary conditions as 1)6. 6|l . By the same computations as in 
Sectional we get that the dynamical large deviation functional is 

Mt]M = I rdt{VH{t),nm+n{t)]VH{t)) (6.8) 



where H has to be obtained from the path tt by using equation 1)6.7(1 with 9(t) 
replaced by n{t). 

This leads to the following Hamilton-Jacobi equation for the quasi potential 

where SV/S9 vanishes at the boundary and 9{±1) — 9±. 
We look for a solution of the form 

6V 9 F 
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By the same computations as in Sectional we reduce the Hamilton-Jacobi equation 
(EHI) to 

, F{1 + F)AF + {9 - F){VFf) =0 (6.11) 



.F^{l + F)^ 

We thus obtain a solution of (|6.9(l considering the functional 

vm ^ ,u {«(.) log ||i + II + «(u)i log Ltm _ } 

where -F(u) has to be computed from 6{u) as the unique strictly increasing solution 
of the boundary value problem 

F" 

F{l + F) ^ + 9~F^0 

{F') (6.13) 

F(±i) = e± 

As for the KMP process it is possible to check that this is the right solution of the 
Hamilton-Jacobi equation H6.9|l . 

By the change of variable F = e*^, it is easy to verify that the right hand side of 
(16.12(1 is strictly convex in (p. We therefore have, analogously to the KMP process, 



V{9) = inf du {9{u) log ^ + [1 + 0{u)] log 
F 7-1 1^ F{u) 



1 + F{u) , F'{u) 
log- 



l + 9iu) ^[9+~9^]/2- 



where the infimum is carried out over all strictly increasing functions F satisfying 
the boundary condition F{±1) — 9±. 



7. Conclusions: few comments on generic models 

For the SEP, the derivation of the rate function for the stationary non equilibrium 
state obtained in [10,11] depends heavily on the details of the microscopic process. 
On the other hand, the variational approach in [5] depends only on the macroscopic 
transport coefficients, bulk diffusion D and mobility a of the system. These are not 
independent functions, they are related by the Einstein relation D{p) = g{p)x{p)^^ 
where x(p) is the compressibility. It is defined as x{p)~^ = ^'{p) = /o (p) where 
/o is the (equilibrium) Helmholtz free energy of the system and A is the chemical 
potential. This means in particular that while the derivation in [10,11] is only valid 
for nearest neighbor jumps, the result holds for the general SEP. In this paper we 
have discussed a model, the KMP process (in fact two models if we consider also its 
dual process), in which the rate functional has an expression very similar to the one 
for the SEP. Here we discuss what are the essential features of the functional form 
of these coefficients in the derivation of the rate functional S. In this discussion we 
shall consider D and a as given and discuss the large deviations properties of the 
nonequilibrium state. 

We discuss only one-dimensional (symmetric) diffusive system with a single con- 
servation law and particle reservoirs at the boundary. Here it will be convenient 
to think of the conserved quantity as the density of particles. For general models, 
the hydrodynamic equation is expected (proven for equilibrium models in [26] and 
in [14, 15] for nonequilibrium under the so called gradient condition)^ to be given 
by a nonlinear diffusion equation with Dirichlet data at boundary, i.e. 

f dtp{t,u) = ]^v{D{p{t,u))Vp{t,u)) 

j p(i,±l) = p± (7-1) 
I P(0,'«) = po{u) 
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where the bulk diffusion D{p) — a{p)x{p) ^ is given by a Green-Kubo formula, 
see e.g. [25, II. 2. 2]. For the KMP process, as well as for the SEP, we simply have 
D = l, i.e. a = x- 

The probability of a large deviations from the hydrodynamic behavior are ex- 
pected (to our knowledge for open systems this has been proven only for the SEP 
in [6], see however [25, II. 3. 7] for an heuristic derivation for equilibrium lattice gas 
models) to have the form (|3.16|) - (|3.17(l where the dynamical cost J[o,t] should be 
of the form 



1 



T 



J[o,T]i7r) ^ - I dt{VH{t),a{^{t))VH{t)) (7.2) 







in which the perturbation H has to be chosen so that the fluctuation tt solves the 
perturbed hydrodynamic 

1. 




V L»(7r(t))V7r(t) - V f7(7r(t)) Vff (t) 



(7.3) 



and u{ti) is the mobility of the system. For the SEP process we have (T(7r) — 7r(l— tt) 
(note that in this case we have < tt < 1) while for the KMP process, respectively 
its dual, we have <j{Tr) — tt^, respectively a{Tr) = 7r(l + tt). 

We first mention the few examples in which it is possible to obtain the rate 
function 5 in a closed form. The following models are however even simpler than 
the SEP or the KMP process since they do not exhibit the non-locality of S, which 
reflects, at the large deviation level, the long range correlations of the system which 
are expected [1,24] to be a generic feature of nonequilibrium models. 

The easiest example is provided by independent particles. In this case we have 
D constant and a linear. The nonequilibrium state is a product measure and it is 
easy to verify that S{p) — j\du f (^p{u) , piuj) , where 

/(p, r) = /o(p) - /o(r) + (p - r)/^(r) = p\og{p/T) - (p ~ r) (7.4) 

and p is the stationary solution of 1)7. l|l . Another example is the so called zero range 
process, see e.g. [25, II. 7.1]. In this case D{p) = ^'(p) and a{p) = $(/o) where the 
(increasing) function $ depends on the microscopic rates. As shown in [8] the 
nonequilibrium state is again a product measure and, as discussed in [4,5], its 
rate function is S{p) = j^^du f [p{u) , p{u)) for / again given by (|7.4|l with the 
appropriate /q. . These examples (the first being a special case of the second) are 
characterized by the fact that a[p) = C exp{\{p)} where C > is a constant and A 
is the chemical potential. The Einstein formula then gives D{p) — <j'{p). The last 
example is the Ginzburg-Landau model, see e.g. [25, II. 7. 3], where a is a constant 
while D is determined by the Einstein relation. In this case the nonequilibrium 
state is still a product measure and its rate function has the same expression as in 
the zero range process. 

We note that for the SEP, as well as for the KMP process and its dual, we have 
D{p) constant and ^{p) a second order polynomial in p. We next show that an 
expression of the rate function S of the nonequilibrium state can be derived under 
a general hypothesis. More precisely, we shall assume that the diffusion coefficient 
D{p) and the mobility a{p) satisfy the following condition. There exists a constant 
a G M such that for any p ^ t 

a{p) - cr(r) a'ir) 



J^'drDir) D{t) 



+ a{p-T) (7.5) 
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This condition, of course, identifies a rather tiny class of models that in fact coin- 
cides with the class of all the examples discussed. We should not expect to be able 
to obtain S in almost a closed form for any model. As we shall see, the locality of 
the functional S corresponds to the special case (in this class) a — 0. 

Let us first discuss which functions D and a satisfy condition H7.5[l . We rewrite 
it with p and t exchanged 

cr(T) - a{p) _ a'{p) 



/; dr D(r) D{p) 
This equation together with l|7.5|l imply 

a'(r) a'{p) 



+ a{T - p) 



D{t) D{p) 

It is easy to see that this is equivalent to 

a'{r) 



2a(r - p) 



= 2aT + c (7.6) 
D{r) 

with c an arbitrary constant. Condition H7.6|l is a necessary condition for the 
validity of (|7.5|l . We rewrite H7.6|l in the integrated form 

dr rD(r) + c dr D{r) 

and substitute it inside (|7.5|l . We thus obtain 

2a ffdr rD{r) , ^ 

A pair (cr, D) is a solution of H7.5|l if and only if is a solution of (j7.6l) and H7.7(l . 

When a = 0, equation (|7.7|l is always satisfied and 1)7. 6|l becomes cr'(r) = cD{r). 
If c ^ we have the solutions corresponding to zero range dynamics (with an extra 
multiplicative factor c); if c = we have the solutions corresponding to Ginzburg- 
Landau models. 

When a 7^ 0, equation (|7.7() becomes 

p rp 
dr rD{r) = {p + t) / dr D{r) (7.8) 

We differenciate with respect to p and obtain 

(p-t)D{p) = / drD{r) 



that is satisfied if and only if D is constant. Now condition 1)7.6(1 imposes that cr(p) 
is a second order polynomial in p. In this class of solutions fall the simple exclusion 
model, the KMP model and its dual. 

To write the rate functional S we need to introduce a little more notation. We 
let d{p) = fodr D{r), since D > the function d is strictly increasing; we denote 
its inverse by b. We finally set A{ip) :— a{b{(p)) . We next denote partial derivatives 
by a subscript. Let us introduce a function of two variables / = /(p, t) such that 
/pp(p, r) = D{p)a{p)~^ = x~^{p) and normahze / so that f{-,T) has a minimum 
at r and f{T,T) — 0. Therefore 



p 



/(p,t) ^ dr dr' 



1 



X{r') 

It is easy to verify that in the equilibrium case, p-\- — p- ~ pa^ the rate function 5*0 
is simply given by S'o(p) = J^^du f (^p{u) , po) ■ To obtain the rate function S in the 
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nonequilibrium case p+ ^ p„ we introduce the functional of two variables 

Qip, ,) :^ fju [muimu))) - i log ^g)^ } (7.9) 

where p is the equilibrium profile. Note that Vd{p{u)) — D{p{u))Vu is a constant 
since its divergence must vanish in the stationary state. 

We claim that, under condition (|7.5|l . the rate function S can be expressed as 
S{p) ~ Q{p, v[p]) where, given p, the auxiliary function ip — Lp[p\ is the solution of 
the Euler-Lagrange equation SQ{p, ip)/6ip = 0, that is 

1 Ay) _^ p-b{ip) ^ ^ 
d(^(±l))=d(p±) 

where we used /t(/5, t) = —D{t){p — t)/(t{t), b'{ip) = D{b{ip))~^ and the definition 
of A. 

The definition of the functional G and the above equation are not really mean- 
ingful if a = 0, as it is the case for the simple models in which S is local discussed 
above. However, in such a case we understand (|7.1()|l as Aip — whose solution 
is ^p{u) = d{p{u)). Plugging it into the functional G we get, by understanding 
(logl)/0 = 0, the correct local functional S{p). On the other hand, as soon as 
a 7^ 0, the functional 5* is nonlocal. 

To establish the claim, we next show that the functional S solves the Hamilton- 
Jacobi equation 

,SS_ ^(^^^^\ ^ /SS_ 
Sp' Sp / \ Sp 

The argument to conclude the identification of S with the quasi potential, as defined 
in (|4.2|) , is indeed essentially the one carried out in section 14.21 and it is therefore 
omitted. 

By the definition of S, we get SS{p)/Sp = SG{p, v)/Sp — fpip, Kf)) so that the 
left hand side of (|7.11() . after an integration by parts, reduces to 



a{p)f,,(p, b{^))Vp + a{p)^^^^^^V^ ~ D{p)Vp) 

where we used the Einstein relation (T{p)fpp{p,b{ip)) ~ D{p) in the first step and 
fpT{p,b{ip)^ — —D{b{(p))/A{ip) in the second one. 

We next write Vd{p) = V[d(p) — (p] + Vp and integrate by parts the first term 
in the last expression above (recall that d{p) and Lp satisfy the same boundary 
conditions). We finally get that the left hand side of H7.11|l equals 

We therefore find that the functional S solves the Hamilton-Jacobi equation 1)7.11(1 
provided p satisfies the equation 

A{p)Ap + \ A\p) + {VpY = (7.12) 

L d{p) -ip 
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In general, we have no reason to expect to be able to express the solution of the 
functional derivative equation H7.11|l by a boundary value problem analogous to 
(|7.12(l . it simply works under our special assumption. 

Up to this point we did not yet really use condition l|7.5(l but, to complete the 
argument, we need to show that (|7.12|) is equivalent to the Euler-Lagrange equation 
1)7. 10|) . By writing H7.5|l with r = b{ip) we get 

d{p) -ip 

and, by comparing (|7.1()(l with 17.12(1 . we see that they are indeed equivalent under 
the above condition. 

As we emphasized, the rate function for SEP is obtained by taking the supremum 
over (fi of G{p, ip), while for the the KMP process we need to take the infimum. We 
can now realize that this depends on the sign of a. Indeed for a > (as in the 
KMP process) the functional Q {p, (p) is concave in V(p while it is convex for a < 
(as in the SEP). 

It is quite tempting to extend the previous derivation to a broader class of models, 
possibly by a different definition of the trial functional Q, however our attempts in 
this direction were not successful. 
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